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Probing the time-variation of the fine-structure constant: 
Results based on Si IV doublets from a UVES sample * 
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Abstract. We report a new constraint on the variation of the fine-structure constant based on the analysis of 15 
Si IV doublets selected from a ESO-UVES sample. We find Aa/a = (-1-0.15 ± 0.43) x 10~^ over a redshift range 
of 1.59 < 2: < 2.92 which is consistent with no variation in a. This result represents a factor of three improvement 
on the constraint on Aa/a based on Si iv doublets compared to the published results in the literature. The alkali 
doublet method used here avoids the implicit assumptions used in the many-multiplet method that chemical and 
ionization inhomogeneities are negligible and isotopic abundances are close to the terrestrial value. 
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1. Introduction 

Some of the modern theories of fundamental physics, in 
particular SUSY GUT and Super-string theory, motivate 
experimental searches of possible variations in the fine- 
structure constant. Such theories require the existence of 
extra 'compactified' spatial dimensions and allow for the 
cosmological evolution of their scale size. As a result, these 
theories naturally predict the cosmological variation of 
fundamental constants in a 4-dimensional sub-space (Uzan 
2003 and reference therein). 

In the framework of the standard Big-bang model, 
quasar spectra can be used as an important tool to test 
the variation of the fine-structure constant, a — e'^/hc, 
by allowing one to measure its value at different redshifts. 
Bahcall, Sargent & Schmidt (1967) were the first to use 
the absorption lines of alkali-doublets seen in QSO spectra 
to constrain the variation of this quantity. Their analy- 
sis provided Aa/a = (a^ — ao)/ao = (—2 ± 5) x 10~^ 
at a redshift z ^ 1.95. Here ao refers to the value of 
the fine-structure constant on Earth and a, to its value 



at redshift z. Since then several authors have used the 
alkali-doublet method (AD method) to constrain the vari- 
ation of a (Wolfe, Brown & Roberts 1976; Levshakov 1994; 
Potekhin & Varshalovich 1994; Cowie & Songaila, 1995; 
Varshalovich, Panchuk & Ivanchik, 1996; Varshalovich, 
Potekhin & Ivanchik 2000; Martinez, Vladilo & Bonifacio, 
2003) . The method is based on the fact that the separation 
between energy levels caused by fine-structure interactions 
is proportional to a^ with the leading term of energy level 
being proportional to a^. As a result, to a very high ac- 
curacy, the relative separation of a fine-structure doublet, 
(A2 — Ai)/A = AA/A, will be proportional to a^. Here Ai 
and A2 are, respectively, the rest wavelength correspond- 
ing to transition ^Si/2 ^^P3/2 and ^Si/2 ^^Pi/2 of the 
alkali-doublet and A is the average value of Ai and A2. 
Varshalovich, Potekhin & Ivanchik (2000) give the follow- 
ing relation between Aa/a and the values of (AA/A) at 
redshifts and z: 



Aa 

a 



cr 



(AA/A), 



(AA/A)o 



(1) 
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where "cr"(w 1) represents the higher order relativistic 
correction. 

Actually, the dependence of rest wavelengths to the 
variation of a is parameterized using the fitting function 
given by Dzuba et al. (1999a) 



uj = u}Q + qiX + q2y 



(2) 



2 
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Here loq and co are, respectively, the vacuum wave num- 
ber (in units of cm~^) measured in the laboratory and 
in the absorption system at redshift z. x and y are the 
dimensionless numbers defined as x — (a^/ao)^ — 1 and 
y = (az/ao)^ ~ 1- The sensitivity coefficients qi and q2 
are obtained using many-body relativistic calculations (see 
Dzuba et al. 1999a). For a given doublet and Aa/a << 1, 
Murphy et al. (2001) have shown that Eq. (2) reduces to 
Eq. (1) with, 

Sqi + 6q2 

CT ~ 

5qi + 26q2 

The value of "cr" for Si iv is 0.8914 when using the q 
coefhcients as given in Table |21 (see below) . 
The AD method works with emission as well as absorp- 
tion lines. However emission lines are usually broad as 
compared to absorption lines. As a result, the constraints 
obtained from emission lines are not as precise as those 
derived from absorption lines. Bahcall et al. (2004) have 
recently found Aa/a — (0.7 ± 1.4) x 10~^ using O ill 
emission lines from QSOs. The most stringent constraint 
from alkali- doublet absorption lines has been obtained by 
Murphy et al. (2001), Aa/a = (-0.5 ±1.3) x 10"^ by 
analyzing a KECK/HIRES sample of 21 Si iv doublets 
observed along 8 QSO sight lines. 

The generalization of this method known as the many- 
multiplet (MM) method (Dzuba et al. 1999b) makes use 
of a combination of transitions from different species. 
The sensitivity coefficients qi and q2 of heavier ele- 
ments are found to be an order of magnitude higher 
than those of lighter elements. As a result the MM 
method gives an order of magnitude better precision in 
the measurement of Aa/a. Application of MM method 
to KECK/HIRES data resulted in the measurement of 
Aa/a = (-0.57 ±0.10) x 10~^ over the redshift range 
0.2 < z < 3.7 (Murphy et al. 2003). However our re- 
cent investigation (Srianand et al. 2004 and Chand et al. 
2004) using very high quality UVES data and well de- 
fined selection criteria resulted instead in a null detection 
of Aa/a (Aa/a = (-0.06 ±0.06) x 10"^) over the red- 
shift range 0.4 < z < 2.3. 

However, these results based on the MM method hinge on 
two assumptions: (i) ionization and chemical homogeneity 
and (ii) isotopic abundances of Mg ii close to the terres- 
trial value. Even though these are reasonable assumptions 
one cannot completely rule out systematic biases induced 
by them in the analysis, especially when one is looking 
for very small effects. On the contrary one can completely 
avoid the assumption of homogeneity in the case of the 
AD method because, by construction, the two lines of the 
doublet must have the same profile (see also Quast et al. 
2004) . Also the effect of isotopic shifts is negligible in the 
case of Si iv doublets (see Section 3.4 of Murphy et al. 
2001). Therefore it is important to increase the precision 
of Aa/a measurements based on the AD method. This can 
be achieved by (a) increasing the S/N ratio and spectral 
resolution of the data used; (b) increasing the sample size. 
The S/N ratio of the data used by Murphy et al. (2001) 



Table 1. List of Si iv doublets in our sample 



QSO 


Ze.vrL 


^ahs 


Comments 


HE 1341-1020 


2.135 


1.915 








2.147 


saturated 


Q 0122-380 


2.190 


1.906 








1.969 
1.973 








1.975 




PKS 0237-23 


2.222 


1.597* 




HE 0001-2340 


2.263 


2.183 




Q 0109-3518 




2.045 


contaminated 


HE 2217-2818 


2.414 


1.965 


weak & blend 






2.186 


blend 


Q 0329-385 


2.435 


2.251 


contaminated 


HE 1158-1843 


2.449 


2.266 


blend 


HE 1347-2457 


2.611 


2.329 




Q 0453-423 


2.658 


2.276 


contaminated 






2.502 


blend 


PKS 0329-255 


2.703 


2.328 


unstable 






2.454 








2.455 




Q 0002-422 


2.767 


2.167* 








2.301 


contaminated 






2.464 




HE 0151-4326 


2.789 


2.451 








2.493 




HE 2347-4342 


2.871 


2.735 


contaminated 


HE 0940-1050 


3.084 


2.667 


contaminated 






2.828 








2.830 


unstable 






2.916 


blend 


PKS 2126-158 


3.280 


2.727 


contaminated 






2.768 


contaminated 






2.907 


contaminated 


Q 0420-388 


3.117 


3.087 


broad 



'*' Si IV system below Lya emission 

but included in our sample (see text), 
"contaminated" Si iv doublet with inconsistent 

profiles due to contamination of absorption 

lines from other systems, 
"blend" the majority of the components 

have separations less than the 

individual h values, 
"unstable" Si iv doublet with unstable Voigt 

profile fit. 



is in the range 15-40 per pixel and the spectral resolution 
is R 34000. In this paper our motivation is to improve 
the Aa/a measurements by using the alkali doublets 
detected in our UVES data of higher S /N and resolution. 
Si IV is used instead of C iv because wavelengths are bet- 
ter known for Si iv (Griesman & Kling 2000; Petitjean & 
Aracil 2004a) and q coefficients are larger (see Table 2). 
The organization of the paper is as follows. In Section 
2 we briefly describe our data sample and analysis. The 
importance of selection criteria is discussed in Section 3 
and discussion of individual systems are given in Section 4. 
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The results and overall discussion are presented in Section 
5. 



2. Data sample and analysis 

2.1. Data sample 

The data used in this study have been obtained with 
the Ultra-violet and Visible Echelle Spectrograph (UVES) 
mounted on the ESO Kuyen 8.2 m telescope at the Paranal 
observatory for the ESO-VLT Large Programme "Cosmic 
evolution of the intergalactic medium" (PI Jacqueline 
Bergeron). This data set corresponds to a homogeneous 
sample of 18 QSO lines of sight with a Lyman-a rcdshift 
range of 1.7 to 3.2. The detailed quantitative description 
of data calibration are presented in Aracil et al. (2003) and 
Chand et al. (2004, here after Paper I). Briefly, the data is 
reduced using the UVES pipeline. Addition of individual 
exposures is performed by a sliding window and weight- 
ing the signal by the errors in each pixel. The final S/N 
ratio is about 60-80 per pixel and the median resolution 
R ^ 45000. The continuum is fitted using an automated 
continuum fitting procedure (Aracil et al. 2003). 

The Si IV systems detected in our data set are listed 
in Tabled There are 31 Si iv systems redshifted beyond 
the Lya emission line from the quasar. In addition, two 
systems, marked with an asterisk (*) in Table 1, falling in 
the Lyman-a forest have well defined narrow components 
and have therefore been incorporated in the sample. 

Among the 31 systems (beyond the Lya emission), 9 
systems are not considered in the analysis because they 
are contaminated by other metal lines and/or atmospheric 
absorption. They are noted as "contaminated" in column 
4 of Table ^ In addition, 4 other systems are rejected 
for the following reasons. One is completely saturated 
( "saturated" in Table and one is a very broad system 
( "broad" ) . The profile of this broad system is spread over 
350 km s^^, and has a few bad pixels in the central part of 
the Si IVA1402 line. The other two systems (marked with 
"unstable" in column 4 of Table ^ were rejected during 
the analysis as we found that the component structure of 
their best-fit is not stable (as discussed in Section 3.2). 
Thus from a total of 33 systems, we exclude 13 systems (9 
contaminated, one saturated, one very broad and 2 sys- 
tems having large uncertainties in the component struc- 
ture) and are left with 20 systems. 

We have shown in Paper I (using detailed simulations) 
that it is better to avoid weak, or heavily saturated, or 
strongly blended absorption lines in order to obtain a bet- 
ter accuracy on Aa/a measurements. Indeed, such sys- 
tems can induce false alarm detection of non-zero Aa/a . 
For a typical S/N ratio of 70 and a median h Doppler 
parameter of 9 km s~^ as seen in our sample, we define a 
lower limit for A^(Si iv) of 1.55 x 10^^ cm^^ so that both 
lines of the doublets are detected at more than 5tT level. 
As in Paper I, we define a multi-component system to be 
unblended if the majority of its components have separa- 
tions larger than the individual h values. We apply these 
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Fig. 1. Left panels show on a velocity scale Si iv doublet 
data points with error-bars together with the best Voigt- 
profile fit for Aa/a = over plotted as a solid curve. 
The dotted vertical lines mark the position of compo- 
nents. Right panels show the variation of as a func- 
tion of Aa/a for the systems in the corresponding left 
panel. Dark rectangles with error bar indicate the mea- 
sured values of Aa/a with one sigma error-bar obtained 
using Xmin + 1 statistics. Name of QSOs and Zabs are 
stated explicitly. 

Table 2. Atomic data of Si iv doublet used in our 
analysis: 



Ion 


A"(A) 


uio{cm ^) 




92 f° 








(cm^l) 


(cm""-) 


Si IV 


1393.76018(4) 


71748.355(2) 


766 


48 0.5140 


Si IV 


1402.77291(4) 


71287.376(2) 


362 


-8 0.2553 



(a) Griesmann U. & Kling R. (2000) 

(b) Dzuba et al. (1999a) 

(c) Martin & Zalubas (1983) and Kelly (1987) 
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Fig. 2. The first and second panels (from left) show on a velocity scale Si iv doublet data points with error-bars 
together with the best Voigt-profile fit for Aa/a = over plotted as a solid curve. The dotted vertical lines mark the 
position of components. The third and forth panels show respectively the variation of as a function of Aa/a for 
the systems in the first and second panels. Dark rectangles with error bar indicate the measured values of Aa/a with 
the one sigma error-bar obtained using -I- 1 statistics. Name of QSOs and Zabs are also stated. 



criteria after the Voigt profile decomposition of the dou- 
blets (as described in the following Section). Based on the 
best fitted parameters, we decide whether a given system 
satisfies our selection criteria or not. Out of 20 systems for 
which we have performed Voigt profile fitting, 4 systems 
are blended and one system is both blended and weak. We 
mark these system respectively by 'blend' and 'weak' in 
column 4 of Table ^ As a result we are finally left with 
15 Si IV doublets that are useful for Aa/a measure- 
ments. The procedure used for the Aa/a measurement is 
described in the next section. 



2.2. Analysis 

We first carry out a Voigt profile fit for each system assum- 
ing Aa/a = 0. The rest wavelengths for the Si iv doublet 
as well as the other atomic parameter used in the fits are 
summarized in Table[21 The Voigt profile fit is carried out 
by simultaneously varying the column density, N , Doppler 



parameter, 6, and redshift, z, for each component till the 
reduced of the fit is ~ 1. This gives us the required 
number of components to be used in the fit and an initial 
guess oi N , b and z for each component. 

To constrain the variation of a, we use the analytical 
expression of the wave number (cj) as a function of a^/ao 
given by Dzuba et al. (1999a), 



W = Wo + 



1 



92 



1 



(3) 



The sensitivity coefficients qi and q2, laboratory wave 
numbers (wp) and rest wavelengths (A) are listed in 
Table HJ The oscillator strengths (/) used in the analy- 
sis are given in the last column. 

We then consider a change in Aa/a and fit the system 
using the new rest wavelengths as given in Eq.(3). We 
vary Aa/a from -20.0 x lO''^ to 20.0 x 10"^ in steps of 
0.1x10"^ and minimize for each value Aa/a by varying 
N, b, and z. The value of Aa/a at which is minimum 
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z ^ =2.834 toward HE 0940-1050 

abs 

rn ' \ ' \ ' I ' rrx"' \ ^ 




-40 -20 20 40 -8 -4 4 8 

Relative Velocity (km/s) Aa/a(in unit of 10^) 

Fig. 4. The left panel shows the Voigt profile best-fit to the system (solid line) along with sub-components profile 
(dashed lines). The right panel shows the variation of the a-s a function of Aa/a . The curve is median smoothed 
over a smoothing scale of 1 x 10~^. The large fluctuation in the curve show that the minima is local and as a result 
the system is not considered in the final result. Note that the instability of the fit mainly comes from the poorly 
constrained structure of the central feature at ~5 km s^^. 



(Xmin) is accepted as the measured value of the Aa/a from 
the system, provided the reduced of the fit is also ^ 
1. The ncr error is obtain from value of Aa/a where 
= ^ Xmin = assuming the error on data are 
normally distributed (Press et al. 2000, p. 691, Theorem D 
and Appendix A). To be on the conservative side we use as 
1(7 error-bar, the larger of the two values of Aa/a derived 
using — Xmin = ^ from the left and right side of Xmin- 
More details about the validation of our fitting procedure 
using the simulated data set can be found in Paper I. 

In principle we can explicitly use Aa/a as one of the 
fitting parameters (like iV, b and z) in our minimiza- 
tion. However in that case we will have no way of ensuring 
whether the best fitted value of Aa/a is a true value or an 
artifact of local minima due to inconsistent fitting of the 
line profiles. Therefore to ensure that the minima with 
respect to our crucial parameter Aa/a is not a local min- 
imum, we have preferred to use the x^ versus Aa/a curve 
method rather than a single minimization that would si- 
multaneously vary all four parameters, Aa/a , N, b and 
z. On the contrary, we minimize x^ by varying iV, b and z 
for each value of Aa/a . We then determine the minima 
of the x^ versus Aa/a curve. Note that the two methods 
are equivalent as shown by Press et al. (2000, p. 691). For 
the sake of completeness we give details in Appendix A. 

The Voigt profile fits toll individual systems and the 
corresponding plots giving x^ as a function of Aa/a are 
presented in Fig. and The Voigt profile fits to the 
other four systems, resulting in much more precise mea- 
surements of Aa/a , are presented in more detail in Fig.|Sl 
and Fig.|Hl(see Sect. 4). The x^ curves presented for all the 
systems are median smoothed over Aa/a = 0.9 x 10~^, 
to avoid the obvious local fluctuations caused by one or 



two points in the curve. Since the smoothing scale we use 
is about one-half of the limiting accuracy from individual 
Si IV doublets, the flnal result is not affected by the proce- 
dure. The results of the Voigt profile fits are summarized 
in Tableini In this table. Cols. 1 and 2 give the QSOs name 
and the emission rcdshift (zcm)- Cols. 3, 4 and 5 give, re- 
spectively, the mean absorption redshift (zabs) of all the 
components in the system, the measured Aa/a value and 
the reduced x^ of the best-fit. The description regarding 
the component structure is provided in columns 6 to 9. In 
column 6 gives Zabs for individual components, while Cols. 
7 and 8 respectively list the column density and velocity 
dispersion. Column 9 refers to the velocity of individual 
components relative to Zabs (listed in column 3). 

In the course of the analysis we also noticed a few 
systems having weak components on the edge of stronger 
components (a specific example is discussed in Sect. 4.1). 
The parameters (TV, 6, z) of the weak components are 
difficult to constrain from an overall fit. We therefore froze 
alternatively one of the parameters (among N, b and z) of 
the weak component and found that the final Aa/a does 
not depend much on the choice of the frozen parameter. 
To be on the conservative side we have taken as the final 
error the largest error of all determinations. Such systems 
have zero errors for the corresponding frozen parameter in 
Table 3. 

3. Importance of the selection 

3.1. Errors versus strength of the absorption 

As an illustration of the importance of the selection of 
the systems, we plot in Fig. O Aa/a measured from the 
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Fig. 5. Left panels shows on a velocity scale Si iv doublet data points with error-bars together with the best Voigt- 
profile fit for Aa/a = over plotted as a solid line. The Voigt profiles of individual sub-components are shown by 
dashed lines. Right panels show the variation of as a function of Aa/a . The dark rectangle with error bars indicates 
the measured value of Aa/a with one sigma error-bars derived using the criteria A^^ = 1 around the minima. The 
name of QSOs, Zabs and the value of measured Aa/a are also stated explicitly. 



20 Si IV systems in our sample versus the column density 
of the strongest component in the system. Open circles 
are for systems that we define as "blended" (they do not 
pass our selection criteria) and filled circles are for systems 
that pass this selection criteria. It is apparent that errors 
are larger for weaker systems and for blended systems as 
expected from the simulations presented in Paper I. The 
column density of the strongest component is considered 
rather than the total column density since the precision 
in Aa/a measurements is most often dominated by the 
strongest component. 



3.2. Precise measurement versus local minima 

We illustrate here the importance of the selection crite- 
ria with an example of a spurious precise measurement 
which results from local minima caused by an unstable 
fit. This happens in the rejected system at 2;abs= 2.8347 
system toward HE 0940-1050. The Voigt profile best-fit 
(along with the profile of sub-components) and the varia- 
tion of as a function of Aa/a for this system are shown, 
respectively, in the left and right panels of Fig.0] The sys- 
tem is fitted with 5 components and xl — 0.97. The right 
panel of Fig. 0] shows that the curve giving for this 
system possesses a large local fluctuation. The cause of 
such fluctuation is that the fit is unstable, either due to 
the uncertainty in the overall component structure or to 
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Fig. 6. Same as Fig. |S| In addition, the position of sub-components is also marked by tick marks. 



abnormal fluctuation in a few pixels (see Fig. ^ around 
V = n km s^^ in Si ivA1402 ). The uncertainty on the 
measurement as derived from this curve using our proce- 
dure would be underestimated. This kind of measurement 
could eventually dominate the whole statistics of weighted 
mean in the final result. It is important to mention that 
if we had not used the versus Aa/a method, we may 
have considered this local minimum as a true minimum 
(Fig.©. 



4. Notes on individual systems 

Table O shows that some of our measurements have error- 
bars on Aa/a less than or comparable to 2 x 10~^ which 
is much smaller than the average. Such systems need to be 
discussed in more detail because they will dominate the 
weighted mean of the Aa/a measurements. 



4.1. Zabs= 2.464 system toward Q 0002- 422 



This system is spread over a velocity range of about 240 
km s~^. Here we have used only the well detached red 
part of the system, as the profile of the Si ivA 1393 line in 
the blue part is affected by a bad pixel. The best-fit Voigt 
profile and the profiles of its sub-components are shown 
in the top left panel of Fig. |S1 As is evident from the fig- 
ure this system is well fitted with eight sub-components 
{x1 = 1.13). The red-most component of this system 
(V« 52 km s""'^) is found to make the Voigt profile fit 
of the system unstable, if we vary all its parameters at the 
same time. As a result we have frozen one of the param- 
eters to their best-fit value obtained assuming Aa/a —0. 
We performed Aa/a measurements for the three cases of 
freezing one of its parameter among N, b, and z (as dis- 
cussed in Sect. 2.2). We find that all measurements are 
consistent with one another and accept the measurements 
with largest error (i.e Aa/a = (1.00 ± 2.08) x 10"^). The 



8 Hum Chand, Patrick Pctitjean, Raghunathan Srianand & Bastien Aracil: fine-structure constant 

Table 3. Results of the Voigt profile fits to the 15 Si IV systems those satisfy our selection criteria. 



Z(j(,g Aa/a Xu Component Structure 

(mean) (in units of 10~^) Zg-bs logio(N)cm~^ b(km s~^) V(kms~^) 



HE 1341-1020 


2.135 


1.91534 


-6.80±12.54 


0.44 


1 


.915491 


± 


0, 


.000800 


12, 


,447 


±0, 


,023 


8, 


,669 


± 


0, 


,632 


16, 


,530 












1 


.915210 


it 


0, 


,003547 


11 


,799 


±0, 


,145 


5, 


,279 


it 


2, 


,652 


-13, 


,383 












1 


.915319 


± 


0, 


,002501 


11 


,751 


± 0, 


,162 


3, 


333 


± 


2 


,404 


-2 


,171 


Q 0122-380 


2.190 


1.90870 


8.00 ± 7.08 


0.80 


1 


.908656 


± 


0, 


,000725 


12 


,597 


± 0, 


,027 


2, 


,651 


± 


0, 


,561 


-4 


,586 












1 


.908745 


± 


0, 


,000940 


12 


,400 


± 0, 


,042 


2, 


,863 


± 


0, 


,867 


4 


,598 






i.yo9Di 


0.00 It 10. 64 


0.52 


1 


ncn A nc 
.969496 


it 





nnnnco 
,000963 


12 


,364 


it 0. 


,026 


8, 


,185 


± 


0, 


,715 


— 11 


,766 












1 


.969729 


± 


0, 


,003927 


11 


,851 


±0, 


,098 


9, 


,433 


it 


2 


,975 


11 


,766 






1.97337 


-0.30 ± 2.34 


0.87 


1 


.973533 


± 


0, 


,001488 


13, 


,188 


±0, 


,028 


11, 


,869 


it 


0, 


,000 


16, 


,779 












1 


.973729 


± 


0, 


,000358 


12 


,547 


± 0, 


,021 


3, 


,191 


it 


0, 


,326 


36, 


,552 












1 


.973096 


± 


0, 


,001843 


12 


,730 


± 0, 


,050 


7, 


,905 


± 


0, 


,873 


-27, 


,315 












1 


.973275 


± 


0, 


,001371 


13, 


,207 


± 0, 


,186 


4, 


,122 


± 


1 


,780 


-9 


,261 












1 


.973192 


± 


0, 


,001700 


12 


,922 


± 0, 


,114 


2, 


,735 


± 


0, 


,841 


-17, 


,633 












1, 


.973377 


it 


0, 


,002469 


13, 


,543 


±0, 


,083 


8, 


,025 


it 


1 


,193 


1 


.034 






1.97587 


-0.30 ± 10.30 


0.68 


1 


.975917 


± 


0, 


,002038 


12 


,303 


±0, 


,214 


4, 


,522 


± 


1 


,177 


4 


.711 












1 


,975823 


± 





,012910 


12 


,064 


± 0, 


374 


8 


132 


± 




,132 


-4 


.759 


PKS 0237-23 


2.223 


1.59671 


-0.60 ± 2.32 


0.82 


1 


.595925 


± 





,001055 


11 


,763 


± 0, 


()98 


1 


160 


± 


2 


,414 


-90 


.264 












1 


.595902 


± 


0, 


,001960 


12 


,464 


± 0, 


,060 


10, 


,260 


it 


0, 


,851 


-92, 


.922 












1 


.596053 


± 


0, 


,001286 


11 


,703 


±0, 


,112 


3, 


,062 


± 


1 


,549 


-76, 


.473 












1 


.597202 


it 


0, 


,000911 


12, 


,582 


±0, 


,159 


6, 


,297 


it 


1 


,082 


67, 


.265 












1 


.597384 


± 


0, 


,001517 


12 


,736 


± 0, 


136 


8, 


762 


± 


1 


,948 


78, 


.296 












1 


,597151 


± 





,048581 


12 


,253 


± 0, 


704 


21.079 ± 15 


,835 


61, 


.385 












1 


,597293 


± 





,001307 


11 


,873 


± 0, 


256 


1 


650 


± 


2 


,934 


67, 


.788 












1 


,597504 


± 





,003343 


11 


,676 


± 0, 


,357 


3. 


785 


± 


2 


,695 


92 


,165 












1 


,596531 


± 





,000759 


12 


,104 


± 0, 


,041 


4, 


,919 


± 


0, 


,704 


-20 


,245 












1 


,096669 


it 





nnn A nc 
,002426 


11 


,898 


it 0, 


,077 


8, 


,786 


± 


2 


,127 


— 4 


,311 












1 


,596154 


it 


0, 


,000337 


11 


,013 


±0, 


,202 


0, 


,276 


it 


0, 


,128 


-63, 


,807 


HE 0001-2340 


2.263 


2.18469 


-0.80 ± 1.80 


1.14 


2 


,183055 


± 


0, 


,000661 


13, 


,121 


± 0, 


,063 


4, 


,817 


± 


0, 


,352 


-153, 


,733 












2 


,183171 


± 


0, 


,004262 


12 


,855 


± 0, 


,133 


9, 


,143 


± 


2 


,517 


-142 


,796 












2 


,183524 


± 


0, 


,002198 


12 


,453 


± 0, 


,184 


5, 


,042 




2 


,232 


-109, 


,556 












2. 


1 ono 


it 


0. 


nn A nnn 
,004220 


12 


,367 


i 0, 


,224 


6, 


,609 


± 


3, 


,394 


— 121. 


,415 












2, 


,183642 


it 


0, 


,002869 


11 


,853 


±0, 


,231 


2, 


,758 


it 


3, 


,389 


-98, 


,439 












2 


,183769 


± 


0, 


,003596 


11 


,654 


± 0, 


212 


4 


,851 


± 


3, 


,711 


-86, 


,468 












2 


,184555 


± 





,001184 


12 


,534 


± 0, 


()46 


6 


542 


± 


0, 


,987 


-12 


,420 












2 


,184685 


± 


0, 


,001782 


12 


,204 


± 0, 


,084 


4, 


,511 


± 


0, 


,993 


-0 


,180 












2, 


,184380 


it 


0, 


,002295 


11 


,676 


±0, 


,147 


3, 


,821 


± 


2, 


,283 


-28, 


,905 












2, 


,186058 


it 


0, 


,000353 


12, 


,502 


±0, 


,012 


5, 


,220 


it 


0, 


,268 


129, 


,163 












2 


1 r A n 


i 





nn 1 nn A 
,001094 


12 


,021 


i 0. 


035 


6 


257 


± 





,893 


108 


,613 












2 


,187051 


± 





,000462 


13 


,158 


± 0, 


008 


12 


018 


± 





,230 


222 


,706 












2 


,186770 


± 





,003257 


12 


,290 


± 0, 


()62 


13 


227 


± 


2 


,068 


196, 


,226 












2 


,186462 


± 





,001357 


11 


,922 


± 0, 


044 


5. 


964 


± 


0, 


,965 


167, 


,209 












2 


,183988 


± 


0, 


,002104 


11 


,575 


± 0, 


219 


0, 


739 


± 


1 


,589 


-65, 


,828 












2 


,185035 


± 


0, 


,002420 


12 


,100 


± 0, 


,052 


11, 


,632 


± 


1 


,599 


32, 


,790 












2, 


,185417 


± 


0, 


,003240 


12, 


,276 


±0, 


,044 


18, 


,832 


± 


2, 


,403 


68, 


,770 












2 


,184102 


it 


0, 


,011745 


11 


,914 


±0, 


,148 


20, 


,314 


it 


7, 


,976 


-66, 


,092 


HE 1347-2457 


2.611 


2.32918 


-5.60 ± 7.84 


0.98 


2 


,328624 


± 


0, 


,001374 


12 


,286 


± 0, 


,029 


10, 


056 


± 





,985 


-50, 


,424 












2 


,328948 


± 


0, 


,000556 


12 


,616 


± 0, 


,014 


8, 


,712 


± 


0, 


,410 


-21 


,227 












2, 


,329217 


± 


0, 


,001300 


12, 


,111 


±0, 


,040 


7, 


,046 


± 


0, 


,970 


3 


.008 












2, 


,329599 


it 


0, 


,003715 


12, 


,338 


±0, 


,058 


16, 


,509 


it 


1 


,703 


37, 


.426 












2 


,329531 


± 


0, 


,000592 


12 


,281 


± 0, 


,050 


3, 


,913 


± 


0, 


,611 


31, 


.303 


PKS 0329-255 


2.685 


2.45465 


-2.00 ± 4.07 


0.47 


2 


,454707 


± 





,001053 


12 


,579 


± 0, 


149 


3, 


141 


± 


1 


,277 


4, 


.886 












2 


,454874 


± 





,001173 


12 


,542 


± 0, 


062 


4 


178 


± 





,598 


19, 


.400 












2 


,454643 


± 





,002560 


13, 


,215 


± 0, 


056 


11, 


463 


± 





,896 


-0, 


.663 












2 


,454378 


± 


0, 


,004649 


12 


,868 


± 0, 


,065 


16, 


,186 


± 


1 


,530 


-23, 


.686 






2.45573 


0.20 ± 5.03 


0.36 


2, 


,455683 


± 


0, 


,001245 


12, 


,593 


±0, 


,177 


3, 


,340 


± 


0, 


,913 


-3 


.706 












2, 


,455768 


it 


0, 


,007809 


12, 


,496 


±0, 


,236 


7, 


,026 


it 


2 


,405 


3 


.684 


t<j UUU2 — 422 


2.760 


2.10010 


U.OO It 2.89 


1.33 


2 


, 16 ^ 1 14: 


i 





,00092a 


12 


,509 


it 


028 


6, 


464 


± 


0, 


,000 


—36 


.100 












2 


,167906 


± 


0, 


,000245 


13, 


,265 


± 0, 


,007 


6, 


,198 


± 


0, 


,128 


-23, 


.592 












2, 


,168234 


± 


0, 


,000308 


13, 


,430 


±0, 


,006 


8, 


,496 


± 


0, 


,164 


7, 


.473 












2 


,168391 


it 


0, 


,000357 


12, 


,985 


±0, 


,026 


4, 


,172 


it 


0, 


,327 


22, 


.328 












2 


,168473 


± 





,001424 


12 


,750 


± 0, 


039 


9 


784 


± 


0, 


,000 


30, 


,094 






2.46401 


1.00 it 2.08 


1.13 


2 


,463220 


± 





,001224 


12 


,068 


± 0, 


,031 


8 


850 


± 





,783 


-68 


,511 












2 


,463558 


it 


0, 


,000494 


11 


,898 


± 0, 


,037 


1, 


,515 


± 


0, 


,791 


-39, 


,232 












2, 


,464161 


± 


0, 


,001020 


13, 


,713 


±0, 


,027 


10, 


,864 


it 


0, 


,673 


12, 


,988 












2, 


,464356 


it 


0, 


,001321 


13, 


,382 


±0, 


,092 


6, 


,834 


it 


1 


,105 


29, 


,878 












2 


,464470 


± 


0, 


,002380 


12 


,977 


± 0, 


142 


5, 


,334 


± 


0, 


,553 


39 


,748 












2 


,463805 


± 





,016700 


12 


,848 


± 0, 


160 


22, 


,392 


± 


4, 


,165 


-17, 


,840 












2 


,463918 


± 





,000331 


13, 


,651 


± 0, 


016 


6, 


757 


± 


0, 


,194 


-8 


,053 












2 


,464605 


± 


0, 


,000000 


11 


,885 


± 0, 


,049 


7, 


,291 


± 


1 


,068 


51 


,455 


HE 0151-4326 


2.740 


2.45140 


1.00 ± 1.92 


0.54 


2 


,451609 


± 


0, 


,000305 


13, 


,190 


± 0, 


,026 


5, 


,991 


± 


0, 


,231 


17, 


,802 












2, 


,451421 


± 


0, 


,001505 


12, 


,585 


±0, 


,113 


10, 


,638 


± 


3 


,157 


1 


.471 












2 


,451696 


it 


0, 


,001122 


12 


,797 


±0, 


,039 


9, 


,311 


it 


0, 


,000 


26, 


.366 












2 


,451250 


± 





,000908 


12 


,519 


± 0, 


074 


5, 


,674 


± 


0, 


,494 


-13, 


,387 












2 


,451043 


± 





,000683 


12 


,228 


± 0, 


017 


7, 


,899 


± 





,444 


-31 


,396 






2.49265 


4.70 ± 5.79 


0.60 


2 


,492705 


± 


0, 


,000453 


12 


,297 


± 0, 


,015 


4, 


,820 


± 


0, 


,297 


4 


,894 












2, 


,492591 


± 


0, 


,002372 


11 


,400 


±0, 


,093 


3, 


,261 


± 


0, 


,000 


-4, 


.916 


HE 0940-1050 


3.084 


2.82831 


-0.30 it 5.46 


0.68 


2 


,827964 


± 


0, 


,003259 


12 


,413 


± 0, 


,087 


7, 


757 


± 


1 


,541 


-26, 


.941 












2 


,828107 


± 


0, 


,000997 


12 


,501 


± 0, 


070 


4 


463 


± 





,618 


-15, 


,731 












2 


,828288 


± 


0, 


,002861 


12 


,493 


± 0, 


154 


4 


490 


± 


1 


,076 


-1 


,551 












2 


,828393 


± 


0, 


,006348 


12 


,267 


± 0, 


235 


5 


232 


± 


1 


,842 


6 


,670 












2 


,828776 


± 


0, 


,005376 


12 


,447 


± 0, 


068 


15. 


967 


± 


1 


,289 


36 


,675 












2 


,828899 


± 


0, 


,000903 


12 


,369 


± 0, 


,078 


6, 


,716 


± 


0, 


,738 


46, 


,316 












2, 


,827730 


it 


0, 


,001076 


11 


,990 


±0, 


,040 


4, 


,449 


it 


0, 


,712 


-46, 


,288 
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7.94 



Look back time (Gyr) 
9.54 10.54 



1 1 .21 



11.68 



"5' 



i3 



20 



Aa/a= ( 0.06±0.06)x lO 5 
based on 23 Mg Il/Fe II 
systems observed with UVES 



-20 




Acx/cx= (0.15±0.43)x lO 5 
based on 1 f5 Si IV systems 
observed with UVES 




o 



o 



o 



o 



^5 



c; 



o 



o 



o 



o 



0.5 



1.5 2 

Kedshif t 



2.5 



3 



Fig. 7. Measured values of Aa/a versus the absorption redshift of the systems. The squares show our determinations 
from the MM method using Mg ii/ Fe ii systems and the narrow shaded region represents the 3a ahowed range. The 
fiUed circles are our present measurements using Si iv doublets from UVES and open circles are the measurements from 
KECK/HIRES data by Murphy et al. (2001). The weighted mean from our 15 systems is Aa/a = (0.15 ± 0.43) x 10"^. 
The 3(T allowed range (—1.14 x 10~^ < Aa/a < 1.44 x 10^^) is shown by the wider shaded region. 



relatively stronger constraint is mainly due to the presence 
of 3 well separated strong components together with good 
S/N ratios (^ 68 per pixel in the nearby continuum). 



4.2. Zabs= 2.451 system toward HE 0151-4326 



The best-fit Voigt profiles and variation of as a function 
of Aa/a for this system are shown in the bottom panels of 
Fig.|31 The system is very well fitted by 5 components with 
xt = 0.54. The measured Aa/a value from this system is 
(1.00 ± 1.92) x 10~^. The good accuracy is mainly due to 
3 well separated strong components and very good S/N 
{'^ 89 per pixel in the nearby continuum). 



4.3. Zabs= 2.1839 system toward HE 0001-2340 



The top panels of Fig. show the best-fit Voigt pro- 
files and variation of as a function of Aa/a for this 
system. The absorption profiles produced by this system 
are spread over 400 km s^^, but most of the components 
are very well separated. The system is fitted by 18 sub- 
components (indicated by tick-marks) with xt = 1.14. 
We have excluded from the fit the velocity range 32 to 
68 km in the Si ivA1402 profile because of the pres- 
ence of several spurious pixels. The measured Aa/a value 
from this system is (-0.80 ± 1.80) x 10"^ The S/N ra- 
tio of the spectrum in the vicinity of this system is about 
52. The presence of a large number of unblended compo- 
nents many of which are narrow and strong increases the 
precision of the Aa/a measurement. 
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Table 4. Summary of results from UVES and HIRES samples: 



Sample 


Numbers 


z 


Aa/a 


(10-^) 




Xw 




of systems 


range 


weighted mean 


mean 


o 




UVES 


15 


1.59-2.82 


+0.15 ±0.43 


f0.57± 1.05 


4.06 


0.29 


HIRES 


21 


2.01-3.02 


-0.52 ±1.22 


-0.12 ± 1.48 


6.80 


0.95 


UVES+HIRES 


36 


1.59-3.02 


-0.04 ±0.56 


f0.17±0.96 


5.76 


0.67 



20 



10 



3 

< 




-10 



12.5 13 13.5 

Fig. 3. The figure shows /S.a/a measured from the 20 
Si IV systems in our sample versus the column density 
of the strongest component in the system. Open circles 
are for systems that we define as "blended" (they do not 
pass our selection criteria) and filled circles are for systems 
that pass this selection criteria. It is apparent that errors 
are larger for weaker systems and for blended systems as 
expected from the simulations presented in Paper I. The 
vertical dashed lines refer to the lower and upper cut off in 
column density to avoid very weak and heavily saturated 
systems. 



4-4. Zabs= 1-59671 system toward PKS 0237-23 

This system falls in the Lyman-a forest but it has an un- 
blended structure with well separated components; as a 
result we have included it in our analysis. The absorp- 
tion profile of Si IV in this system is spread over about 
240 km s~^. The best- fit Voigt profile along with profile 
of the different sub-components are shown in the bottom 
left panel of Fig. The system is fitted with 11 compo- 
nents with xt = 0.82. The average S/N ratio per pixel 
in the neighboring continuum is about 57. The measured 
Aa/a value from this system is (—0.60 ± 2.32) x 10~^. 
The relatively good precision is mainly due to the pres- 
ence of a large number of well separated and strong enough 
sub-components. 



5. Results and discussion 

The detailed description of our individual measurements 
from 15 Si iv doublets is given in Table |3| The summary 
of the results obtained from different samples is presented 
in Table ^ In column 1 "UVES" refers to our sample 
and "HIRES" refers to the KECK/HIRES sample from 
Murphy et al. (2001). Cols. 2 and 3 list, respectively, the 
number of systems in the sample and their redshift cov- 
erage. The weighted mean and mean value of Aa/a are 
listed respectively in Cols. 4 and 5. Cols. 6, 7 give the stan- 
dard deviation of measurements around the mean and the 
reduced around the weighted mean. 
The weighted mean of the measurements is obtained by 
assigning weights (wj) as 1/error^ and the error on the 
weighted mean is computed by the standard equation. 



Error in x„ 



(4) 



Here Xw I'efers to the reduced x'^ of variable Xi around 
their weighted mean The Xw term takes into account 
the effect of scatter in measurements while computing the 
error on the weighted mean Aa/a . The error on the sim- 
ple mean Aa/a ( column 5 of Table 2J) is computed by 
the central limit theorem (cr/y/N ), assuming the individ- 
ual measurements are Gaussian distributed around their 
mean. 

The distribution of our 15 measurements together with 
the 21 measurements of Murphy et al. (2001) is plotted as 
a function of Zabs and look-back time in Fig. {7\ The look- 
back time corresponding to a given redshift is computed 
in the case of a flat universe with fix = 0.7, flm = 0.3 
and Hq = 68 km s^^Mpc^"'^. The measurements from our 
UVES sample using the MM method on Mg ii systems 
(Paper I) are also shown for comparison. The weighted 
mean value obtained from our analysis over the redshift 
range 1.59 < z< 2.82 is Aa/a = (0.15 ± 0.43) x 10"^ 
The 3cr range (-1.14 x 10"^ < Aa/a < 1.44 x 10"^) 
is shown in Fig. 13 as a shaded region. Our result corre- 
sponds to a factor of three improvement on the constraint 
based on Si iv doublets compared to the previous study by 
Murphy et al. (2001). The increased accuracy in our result 
is mainly due to the better quality of the data (S/N ratio 
of about ~ 70 per pixel and R ~ 45000). Combining our 
sample with KECK/HIRES sample results in a weighted 
mean Aa/a = (— 0.04± 0.56) x 10~^ over a redshift range 
of 1.59 < z < 3.02. The small enhancement in the error 
in the weighted mean Aa/a in the case of the combined 
sample is due to higher Xw- 
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Further improvements at higher redshift can be 
achieved using as MM analysis of multiplets from single 
species such as Ni ii or Fc ii using a well defined high qual- 
ity sample (see e.g. Quast et al. 2004). It is also demon- 
strated that OH and other molecular lines can be used 
to improve limits on the variation of a (see for example 
Chengalur & Kanekar, 2003). In addition other constants 
can be constrained in a similar way. Although it is hard to 
make any quantitative prediction theorists estimate that 
variations in the proton-to-electron mass ratio could be 
larger than that of the fine-structure constant by a factor 
of 10 to 50. It is possible to constrain this constant by 
measuring the wavelengths of radiative transitions pro- 
duced by molecular hydrogen, H2. On-going ESO pro- 
grames have been dedicated to this purpose (Ledoux et 
al. 2003, Ivanchik et al. 2002, Petitjean et al. 2004b). 
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Appendix A: Relation between errors S{Aa/a) 
estimated from the versus 
Aa/a curve and the covariance 
matrix \/Cii 

In general to fit a nonlinear function like a Voigt-profile, 
one would define the merit function 



N 

E 



(% - y{xj : «))^ 



and minimize it to get the best-fit value of the parame- 
ters ttfc where k=l, 2, ... M. The minimization of such a 
nonlinear function involves an iterative process. In a given 
iteration the trial values of parameters are improved till 
one reaches minima. 

Let us suppose that with sufficient accuracy, wc can 
express x^i'^) ^ quadratic function near its minimum. 
Then 



1- 



d - a + -a - D ■ a 



where "d" is a row matrix of M elements and "£>" is a 
M X M matrix (Press et al. 2000, p. 675). 

For such a quadratic function to go from the cur- 
rent parameter Ucur to amin, which is the parameter at 



Xmin^ ^ Newton method of minimization (al- 

most equivalent to a variable metric or Hessian matrix 
method). It states that if amin is the parameter at Xmin 
then Vx'^{o,min) = 0. 

Now let us suppose acur is close enough to amin so that 
to second order we can write 

X i^min) X (^cur) {^min ^ctir) ' d 

~^Ti{^niin ^cur) ' ^ ' {^min ^cur) (-^•^) 



where "d" and "_D" correspond to the first and second 
derivative terms evaluated at Qcur- 

The requirement that Vx^(amm) = in Eq.(A.2) will give 
us 



D.Sa = d 



where 5a = a„ 



(A.3) 



and 



daiUk ' 



duk 



dk = - 



In a more familiar form Eq.(A.3) appears as 

aki ■ Sai = Pk 



(A.4) 



where (3k = ^dk, o-ki = \Dki and summation is assumed 
over the repeated indices. Now the whole eff'ort of mini- 
mization is to make jSk vanish in Eq.(A.4). At this points 
5a = Q and hence acur will be equal to amin- 

Now we apply this general approach to our specific case 
of four parameters Aa/a , iV, b and z. Let us suppose that 
amin = [af",a^*",a^™,a^™] are the best-fit value of 
these parameters achieved by varying all four parameters 
and let Xmini'^min) be the minimum value of . This 
parameter set is not different from the one that is obtained 
from the minima of the x^ versus Aa/a curve. This is 
because the versus Aa/a curve is obtained by standard 
X^ minimization, by varying N, b and z, at every value of 
Aa/a . The reason we have chosen this method is that 
the x^ versus Aa/a curve will allow us to avoid the local 
minima. 

Now let us keep fixed the first element of our parameter 
vector with a value close to a™"", and let us vary the 
other three parameters to achieve the x^ minimization. 
The resulting parameter set is Oc = [al, a%, a§, 04] {a\ has 
been kept fixed) and the minimum value of is XminC'^c)- 

Now substituting a^ur by flc in the more general equa- 
tion Eq.(A.4), we get 



5a- a - 5a = Ax^ 



Xmm('^c) Xmini^min) 



(A.l) a -5a = 13 



(A.5) 
(A.6) 



where 5a = amin — Oc- But here we should keep in mind 
that Oc is also the best fit parameter obtained by x^ min- 
imization, except that its first element was held fixed. We 
know that for the best-fit parameters the corresponding 
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elements of /3 should vanish. Therefore imposing all ele- 
ments of /3, except the first, to be zero, Eq.(A.6) becomes 



• 6a 



givmg 






W 



(A.7) 



Sai/Cn = c 

where we have used the fact that a is the inverse of the 
covariance matrix C. Combining Eq.(A.7) and Eq.(A.5) 
we get 

{Saif = Ax'Cii 

and using the fact that Cn — af, where cti is the error-bar 
on the first parameter, we get the required relation 



6ai = ±\^A)^ai. 



(A.J 



Interestingly 5ai (the difference between two measure- 
ments of first parameter) becomes equal to cti (the er- 
ror obtained using a covariance matrix, while varying all 
the four parameters) when Ax^ = 1. Also when we use 
versus Aa/a curve we obtain the error in Aa/a using 
Ax^ — 1. This means that the estimated errors using both 
the methods are identical. 
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